Spin effects in single-electron transport through carbon nanotube quantum dots 
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We investigate the total spin in an individual single-wall carbon nanotube quantum dot with 
various numbers of electrons in a shell by using the ratio of the saturation currents of the first steps 
of Coulomb staircases for positive and negative biases. The current ratio reflects the total-spin 
transition that is increased or decreased when the dot is connected to strongly asymmetric tunnel 
barriers. Our results indicate that total spin states with and without magnetic fields can be traced 
by this method. 
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The electron-spin configurations in quantum dots are 
important not only for understanding electron-electron 
interactions but also for the realization of spin-based 
quantum computing devices in solid state systems^. Re- 
cently, spin states in carbon nanotube quantum dots 
have been revealed in various systems with the Kondo ef- 
feclj^i^ and the simple shell structurea^'^'^i^. In single- 
wall carbon nanotube (SWCNT) quantum dots, two- 
or four-electron shell structures are predicted, depend- 
ing on the relation between the zero-dimensional (OD) 
level spacing (A) and an energy mismatch (5) in the 
two subbands in the SWCNT^^. The two-electron shell 
structure [5 A/2) originates from the twofold spin de- 
generacy of single-particle levels, and the four-electron 
shell structure [5 <^ A) originates from the predicted 
twofold subband degeneracy in addition to the spin de- 
generacy!^. Because the OD level spacing and charg- 
ing energy [Ec ~ e^/Cs :Cs is the self-capacitance of 
a dot) are significantly larger than the interaction en- 
ergies of the on-site Coulomb energy {5U) and the ex- 
change interaction energy (J), it is possible to maintain 
the shell structures, regardless of the number of electrons 
in the doti^. We have clearly demonstrated the two- 
and four-electron shell structures in the closed SWCNT 
quantum dots and investigated the total-spin transition 
by considering the magnetic field evolution of the single- 
electron transport and by performing excitation spec- 
troscopy measurementa^ii^. 

In this study, we carefully investigate the total spin in 
an individual SWCNT quantum dot by using the ratio of 
the saturation currents for positive and negative biases 
in single-electron transport through the dot. When a 
quantum dot is connected to strongly asymmetric tunnel 
barriers, single-electron tunneling reflects the total-spin 
transition that is increased or decrease d^'^i^^'^^ . We find 
that total-spin transitions in magnetic fields, such as the 
higher-spin states {S = 1), can be traced by this method 
and the experimental results are consistent with the pre- 
vious results^. 

The single-electron transport measurements are car- 
ried out by applying source-drain bias voltage (T4d) be- 
tween left and right barriers. We define the positive po- 



larity of V^d when electrons flow from the left reservoir to 
the right reservoir. Then, we define the saturation cur- 
rents of the first steps of the Coulomb staircases for posi- 
tive and negative VsdS as /+ and /_, respectively. When 
the tunneling rates of the left and right barriers of the 
dot (represented by /l and /r, respectively) are signifi- 
cantly different, the current is only limited by the thicker 
barrier with the smaller tunneling rate. For instance, if 
the left barrier is thicker (II ^ /r), the current of the 
single-electron tunneling is approximately determined by 
Fi}- . In this case, the ratio of saturation current reflects 
the total spin states in quantum dots, with the relation- 
ship /+//_ = (25*^+1 + l)/(2S'Ar + l)^^, where Sn is 
the total-spin state with iV-electrons; S'at+i, the total- 
spin state with {N -\- l)-electrons. Both Sjq and Sjq+i 
are for the ground state. It should be noted that single- 
particle levels in a quantum dot are degenerate when 
this argument is applied. For instance, for the ground 
state transition from Sn = to Sn+i = 1/2, the tun- 
neling rate for positive Vsd is twice that for negative Vsd- 
For positive Vsd, either a spin-up or spin-down electron 
can enter the dot and exit rapidly due to the larger tun- 
neling rate of the right thinner barrier^®. On the other 
hand, for negative Vsd, only one electron (spin- up or spin- 
down) can rapidly enter the dot and exit through the 
left thicker barrier. Therefore, assuming that the tun- 
neling rates are independent of spin states, the current 
ratio is /+//- ~ 2. For the ground state transition from 
Sn = 1/2 to Sn+i = 0, only one electron (spin-up or 
spin-down) can enter the dot and create a spin pair in 
a spin-degenerate single-particle level {S = 0) and one 
of the two electrons exits the dot. Since the current is 
limited by the incoming tunneling rate, the current ra- 
tio becomes I+/I- = 1/2. The experimental results for 
this method have been demonstrated in carbon nanotube 
quantum dots at zero magnetic fleldi^, and in lateral 
quantum dots made of GaAs/AlGaAs two-dimensional 
electron gas systemic. 

SWCNT quantum dots were fabricated by depositing 
metallic (Ti) source and drain contacts on top of an indi- 
vidual SWCNT using 50 keV electron beam lithography. 
In our devices, the entire nanotube between the contacts 
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FIG. 1: (a) Gray scale plot of the difFerential conductance 
dl/dVBd as a function of V^d and at _B = T. The num- 
ber n indicates the number of extra electrons estimated from 
the Goulomb diamond around Vg —0.27 V. is the to- 
tal number of electrons in the dot and is suggested to be an 
even number (see text), (b) Coulomb oscillations in the gate- 
voltage range that corresponds to (a). Inset: Addition energy 
(-Eadd) as a function of n. The alternate change in the addi- 
tion energy is observed, indicating the even- odd effecfei^iSa. 
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behaved as a single quantum dot^SiSi. The distance be- 
tween the contacts was designed to be 300 nm, and a 
heavily p-doped Si substrate was used for the apphca- 
tion of the gate voltage Vg. All the electrical transport 
measurements were carried out in a dilution refrigera- 
tor with a superconducting magnet. In our experiments, 
magnetic fields were applied perpendicular to the tube 
axis. 

Figure 1(a) shows a gray scale plot of the differential 
conductance dl/dVsd as a function of Vgd and at i? = 
T. The white-colored regions are Coulomb blockaded re- 
gions, and the lines outside the diamonds are due to the 
effect of zero-dimensional confined levels. Figure 1(b) 
shows the Coulomb oscillations in the region of Fig. 1(a) 
at a small voltage of 14d = 25 ^V. The electron tem- 
perature, Te = 120 mK (^ 10 fJ,eY), is obtained by the 
theoretical fitting of an individual Coulomb peak (not 
shown here) ^^'^'^ . The electrostatic potential of the dot 
(pdot is linearly related to Vg with the standard single- 
particle modelSi. We also deduce that the conversion 
factor is a = d(j)dot/dVg — 0.15 meV/mV from the fit- 
ting of the Coulomb peaks. The distance between ad- 
jacent peaks (AVg), in fact, oscillates as a function of 
the gate voltage with a period of two electrons. It corre- 
sponds to the addition energy i?add = aAVg. The inset 
in Fig. 1(b) shows i?add as a function of the number of 
extra electrons in the dot. The oscillation is clearly ob- 
served, indicating the even-odd effect of the number of 
electrons in the dot^^. These results are similar to those 
of the typical "closed" quantum-dot behavior of an in- 
dividual SWCNT, which has been reported previously^S. 
For larger diamonds, i?add is estimated to be 8 — 10 meV 
when the number of electrons in the dot is even, and 



FIG. 2: (a) Typical Coulomb diamonds (n = 8 to 10). (b) 
I — Vsd characteristics of the intersection of adjacent Coulomb 
diamonds in (a); A*' ^ (A'' -I- 1) (open circles) and (A'' -f 1) ^ 
(A'' -I- 2) (closed circles), (c) Typical Coulomb oscillations 
(n = to 1) as a function of Vg for Vad from —2 to 2 mV. A 
magnetic field of 7.0 T was applied perpendicular to the tube 
axis. Each peak is shifted for clarity, (d) 7— Vsd characteristics 
along the dashed lines (i) (open circles) and (ii) (closed circles) 
in (c). 

for the smaller diamonds, it is estimated to be ~ 6 meV 
when the number of electrons is odd. S is distributed in 
the range of 2 — 3 meV for different Coulomb diamonds, 
which is estimated from the excitation lines outside the 
Coulomb diamonds in Fig. 1(a). 

We proceed to measure the saturation currents (/+ 
and /_) with and without magnetic fields for the re- 
gion shown in Fig. 1(a). Figure 2(b) shows the typical 
/ — V^d characteristics between N {n = 8)- and (A^ -I- 1)- 
electron states (open circles) and between (A^ + 1)- and 
{N + 2)- electron states (closed circles) in Fig. 2(a) at 
S = T. In this case, the ratio of saturation currents 
(/_)-//_) is obtained as ~2 for A^ ^ (A^ + 1) transition 
and ^0.5 for (A^ -I- 1) ^ {N + 2) transition. In magnetic 
fields, the Zeeman splitting of single-particle levels is ex- 
pected, and new current steps appear for each individual 
single-particle level. The current steps of Coulomb stair- 
cases are very small in low magnetic fields and therefore 
the current ratio is not obtained by measuring the these 
staircases. However, the current steps can be clearly ob- 
served by measuring the Coulomb oscillations for large 
positive and negative biases. Figure 2(c) shows the typi- 
cal Coulomb oscillations (n = to 1) as a function of Vg 
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FIG. 3: (a) — (c)Magnetic field dependence of the saturation 
current ratio (/+//_) from n = 4 to 10 in Fig. 1(a) for the 
transition from even to odd-A'^ state (open circles) and for 
the transition from odd to even-N state (closed circles), (d) 
Magnetic field evolution of Coulomb peaks up to 10 T from 
n = 4 to 10 in Fig. 1(a). Kd = 0.1 mV. The black-colored 
regions show the Coulomb blockaded regime, and the white- 
colored regions show the Coulomb peaks. 



when Vsd is changed from —2 to 2 mV a.t B = 7.0 T. Fig- 
ure 2(d) (closed circles) shows the current values plotted 
along the first step of a Coulomb peak (dashed line (ii)) 
as a function of Vgd- The ratio of the saturation current 
can be easily determined even if the width of the current 
step is small. When the currents are plotted along the 
dashed line (i), which corresponds to the usual / — V^d 
characteristics at a fixed value of Vg, the obtained result 
is shown in Fig. 2(d) (open circles). which is 

obtained as ^1 in this case, has the same values for each 
method. We have carefully investigated the spin states in 
the carbon nanotube quantum dot by performing these 
analyses in the magnetic fields from to 10 T. 

Figures 3(a) — (c) show the magnetic field dependence 
of I+/I- from n = 4 to 10 in Fig. 1(a). Open circles 
indicate the current ratio for the transition from an even 
to an odd number of electrons and closed circles indicate 
that for the transition from an odd to an even number of 
electrons. Up to i? ~ 1.5 T, I+/I- is ~2 for all the even- 
to-odd transitions and ^^0.5 for all the odd-to-even transi- 
tions. These results indicate that the total spin changes 
from 5 = to S* = 1/2 when the number of electrons 
changes from even to odd. On the other hand, the total 
spin changes from 5 = 1/2 to 5* = when the number of 
electrons changes from odd to even. Figure 3(d) shows 
the magnetic field evolution of Coulomb peaks up to 10 
T from n = 4 to 10 in Fig. 1(a). The even-odd effect is 
again clearly demonstrated in evolution of the Coulomb 
peaks as a function of the magnetic field. All the peaks 
shift linearly as the magnetic field is increased, and the 
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FIG. 4: (a) — (d) Magnetic field dependence of I+/I- from 
n = to 4 in Fig. 1(a), for the transition from even to odd-A*" 
state (open circles) and for the transition from odd to even- A'' 
state (closed circles), (d) Magnetic field evolution of Coulomb 
peaks up to 10 T from n = to 4 in Fig. 1(a). Vad =0.1 
mV. This plot is obtained by the same method that used to 
obtain the plot in Fig. 3(d). The dashed line connects the 
kink positions of two Coulomb peaks (n = 1 -i-^ 2 and 2^3). 
This result is simillar to the four-electron shell-filling scheme; 
an internal spin-flip may occur in this line^i^. 



direction of the shift changes alternately as the number 
of electrons is increased by one unit^^. This result is ex- 
plained by the Zeeman splitting of each spin-degenerate 
single-particle level. The Zeeman splitting of the single- 
particle levels gives a g-factor of ^2, and the total spin 
changes between and 1/2 as the number of electrons 
is increasedi^ii^. These results agree with the results of 
the current ratio in Fig. 3(a) — (c) at _B = T. Then, at 
B ~ 1.5 T, the value of I+/I- jumps to ~1.0 and re- 
mains constant in the magnetic field larger than B — 2.0 
T. This effect is explained by the fact that the spin de- 
generacies of the single-particle levels are lifted by the 
Zeeman effect. In this situation, the argument that as- 
sumes the spin-degenerate levels is not applicable. Only 
one electron (spin-up or spin-down) can tunnel through 
the (spin-up or spin-down) Zeeman split level, and the 
tunneling rate is independent of the bias polarity. The 
Zeeman energy at B = 2.0 T is ~230 /icV, which is sig- 
nificantly larger than the thermal energy in our experi- 
ment, where the current steps due to the Zeeman effect 
are clearly observed?^. 

Now, we consider the region of rt = to 4 shown in 
Fig. 1(a); this region shows a slightly more complicated 
behavior in magnetic fields. Figures 4(a) — (d) show the 
magnetic field dependence of the saturation current ratio 
from n = to 4 region. The n = <-> 1 and 3^4 
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FIG. 5: Mechanisims of the n = 1 <-> 2 transition for (a) pos- 
itive V^d and (b) negative Vsd and the n — 2 ^ 3 transion 
for (c) positive Vsd and (d) negative Vsd in magnetic fields in 
the single-particle model. A and B sites represent the ground 
state and the first excited state, respectively. Each single- 
particle levels have undergone Zeeman splitting. The right 
barrier is thinner in our experiment; therefore, tunneling cur- 
rent is limited only by the tunneling rate of the thicker left 
barrier (see text). 



transitions (Figs. 4(a) and (d), respectively) produce the 
same results as those in Fig. 3. 

For the n = 1 ^ 2 transition, /+//- changes from 
~1 to -1.5 at S = 5.5 T in Fig. 4(b), which theoret- 
ically corresponds to 3/2 for Sn = 1/2 and Sn+i = 1- 
In magnetic fields, single-particle levels are lifted and no 
longer degenerate; consequently, the formula of the cur- 
rent ratio may not be applicable. However, in this situ- 
ation, two-electron spin states are degenerate due to the 
level crossings between the ^-spin-down state and the B- 
spin-up state, as shown schematically in Figs. 5(a) and 
(b) , and the formula of the current ratio may be applica- 
ble. Therefore, a higher-spin transition can be traced in 
region (I) in Fig. 4(b). When the degenerate levels be- 
tween the ^-spin-down state and the i?-spin-up state dif- 
fer slightly, assuming that the crossing levels are lifted in 
Figs. 5(a) and (b), two electrons can enter the dot when 
Vsd is positive, and two electrons can exit the dot when 
Kid is negative. Therefore, the value of the current ra- 



tio jumps to ~1.0 again for larger magnetic fields (region 

(II) ), and is consistent with the elementary value of the 
single-particle model. On the other hand, for rt = 2 ^ 3, 
/_!_//_ changes from —1 to '^0.67 at i? 6.0 T in region 

(III) in Fig. 4(c). I+/I- agrees with the theoretical value 
of 2/3 for Sn = 1 and Sn+i = 1/2. At B > 8.0 T, /+//_ 
jumps from —0.67 to —0.33 for n = 2 <-^- 3 in region (IV) 
in Fig. 4(c). In this case, three-electron spin states are no 
longer degenerate due to the Zeeman effect; therefore we 
can consider the simple single-particle model, assuming 
that the crossing levels are lifted in Figs. 5(c) and (d). 
Only one electron (spin-down) can enter the dot when the 
Vsd is positive. On the other hand, three electrons can 
exit the dot when the Vsd is negative; therefore, is 
equal to 1 /3 and is consistent with the elementary value 
of the single-particle model. 

The above results indicate that at i? = 6 T, the total 
spin changes from 0— > 1/2— > 1 1/2— >0 as n in- 
creases. Figure 4(d) shows the magnetic field evolution 
of the Coulomb peaks up to 10 T from n = to 4 in 
Fig. 1(a). This behavior is similar to the four-electron 
shell-filling scheme, as reported previously''. These re- 
sults are consistent with the results of the current ratio 
shown in Fig. 4(a) — (d). We show that the higher-spin 
transitions {S = 1) in magnetic fields can be traced using 
this method. 

In conclusion, we have determined the total spin and 
its magnetic field dependence by using the ratio of 
the saturation currents for positive and negative biases. 
Higher-spin transitions in magnetic fields can be traced 
by this method when the single-particle levels are degen- 
erate; then, the single-particle model cannot be applied 
in the sigle-electron transport such that shown in Fig. 
5. Furthermore, these results indicate that the spin-up 
and spin-down states of a single-particle level that has 
undergone Zeeman splitting have an equivalent tunnel- 
ing rate in the SWCNT quantum dot. The experimental 
results are consistent with those of the total-spin study 
due to the magnetic field evolution of Coulomb peaks, as 
reported earlier. This method can be used to determine 
the total spin in quantum dots. 
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